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Abstract 

We compute non-Gaussianities in A/'-fiation, a string motivated model of assisted inflation with quadratic, 
separable potentials and masses given by the Marcenko-Pastur distribution. After estimating parameters 
characterizing the bi- and trispectrum in the horizon crossing approximation, we focus on the non-linearity 
parameter /jvl, a measure of the bispectrum; we compute its magnitude for narrow and broad spreads of 
masses, including the evolution of modes after horizon crossing. We identify additional contributions due 
to said evolution and show that they are suppressed as long as the fields are evolving slowly. This renders 
A/'-fiation indistinguishable from simple single-field models in this regime. Larger non-Gaussianities are 
expected to arise for fields that start to evolve faster, and we suggest an analytic technique to estimate their 
contribution. However, such fast roll during infiation is not expected in A/'-fiation, leaving (p)re-heating as 
the main additional candidate for generating non-Gaussianities. 
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I. INTRODUCTION 

The generation of a nearly scale invariant spectrum of primordial perturbations, as observed 
in the cosmic microwave background radiation (CMBR) [l|, 13, [3, or large scale structure (LSS) 
surveys , is a decisive indicator of an inflationary epoch in the early universe. During inflation 
Q], an accelerated expansion of the early universe, perturbations are stretched beyond the Hubble 
radius and re-enter the late universe, making them observable to us, e.g. in the CMBR. 

A potential approach to discriminating between inflationary models consist of measuring devia- 
tions from purely Gaussian statistics; Most simple models of inflation predict quasi scale invariant, 
nearly Gaussian, adiabatic perturbations. However, more intricate models such as multi-field infla- 
tionary ones (see e.g. [S] for a review), might deviate from Gaussianity. Henceforth, a consideration 
of higher order correlation functions, such as the bispectrum or trispectrum, could potentially shed 
light into the fundamental physics responsible for generating primordial fluctuations P, [l^ . 
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A rough estimate of non-Gaussianities, as measured by the bispectrum, is the non-hnearity 
parameter /nl, properly defined later in the text. Currently, primordial non-Gaussian features in 
the cosmic microwave background radiation have not been detected. For instance, the observational 
bound on the non-linearity parameter obtained from the WMAP3 data set alone is —54 < /nl < 
114 [l|,0,Q], and future experiments will improve upon it [isl . Q, [isl . 11^ . Hence, theoretical 
predictions for /tvl in well motivated inflationary models should be made before the increasingly 
improved observational data is available. 

All models of inflation predict, to a certain extent, some level of primordial non-Gaussianity. 
^n single-field models, non-Gaussianities are small, usually of the order of the slow roll parameters 
l?! because fluctuations freeze out once their wavelength crosses the Hubble radius. In multi-field 
models, however, the presence of multiple light degrees of freedom perpendicular to the adiabatic 
direction leads to the generation of isocurvature perturbations, which in turn allow the comoving 
curvature perturbation to evolve after horizon crossing (HC). This leads to an additional source 
of potentially detectable non-Gaussianity [18,]. 

One such multi-field modelj^ assisted inflation, which was originally proposed to relax fine 
tuning of potentials, (see e.g. 



IS 



2C 



21 



22l|): it relies on M scalar fields, preferably uncoupled, 



which assist each other in driving an inflationary phase. Even though each individual field may 
not be able to generate an extended period of inflation on its own, they can do so cooperatively. 
This phenomenological model is attractive, since super Planckian initial values of the fields can be 



avoided 



IS 



20 



2l| if the number of fields is large. A string motivated implementation of assisted 



inflation is A^-flation 



251, 



27 



291, 



30 



3ll |: here the many fie^ 



ds are identified with axions 



3l| 



arising from some KKLT compactification of type JIB string theory 

Although the existence of AA-flation via a concrete construction has not been proven, it pro- 
vides a test-bed for multi-field inflation. In [31] it is argued, but not proven, that an effective 
quadratic potential for each field without cross couplings may be attainable. Based on results 
of random matrix theory, it is further argued in [si] that the masses for the M fields, originally 



considered identical in 



23], would conform to the Marcenko-Pastur (MP) distribution 



33]. This 



spectrum of masses is controlled by only two constants: the average mass and a parameter fixing 
the width/shape of the spectrum, which may be identified with the ratio of the number of axions 
to the total dimension of the moduli space in a given construction. The effect of this distribution 



^ See also [23l . [23 ] for another implementation of assisted inflation within M-theory, making use of multiple M5- 
branes. 
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on the power-spectrum has been computed in [25j, |31[, yielding a sHghtly redder spectral index. 



Non-Gaussianities arising in multi-field inflation have been considered recently: the starting 
joint is the three-point correlation function just after horizon crossing, which was first derived in 

d 2. 



34 



33] • Thereafter, one can evolve non-Gaussianities via the (non-linear) 5A^-formalism 
If the horizon-crossing (HC) approximation is used, models of assisted inflation become indistin- 
guishable from their single-fleld analogs This is easily understood since, by discarding the 
evolution of modes after horizon crossing, the main distinguishing feature of multi-fleld models is 
neglected. As a result, the spectrum of scalar and tensor perturbations is exactly the same ^ to 



the one generated by an effective single-field model 

In this paper we focus on AA-flation with an arbitrary, but large number of fields, so that the 
Marcenko-Pastur distribution can be used. We consider narrow and broad mass spectra, extending 
the study of SeJ where the formalism employed in this study was developed; its application, how- 
ever, was limited to simple toy models. Here, we analytically compute the non-linearity parameter 
Inl without assuming the horizon crossing approximation, but within the slow roll approximation. 
A comparison with simple estimates as well as the HC-limit, which is re-derived for completeness, 
reveals that additional contributions due to the evolution of modes after horizon crossing are 
present, but their magnitude is limited to a few percent of the HC result. Hereafter, we relax the 
slow roll approximation and argue that large, but possibly transient contributions to f^L should be 
expected from faster rolling fields, which are however not expected during AA-fiation, since heavy 
fields should evolve slowly up until (p)re-heating commences. Nevertheless, an analytic method 
to retain the main physical effect of fields starting to roll faster during inflation is suggested: an 
effective single-field model with steps in its potential. 

The concrete outline is as follows: in Section |TI] we review AA-fiation, focusing on the evolution 
during slow roll. After that, we introduce non-Gaussianities in Section IIIII and provide, for later 
reference, the horizon crossing results of the non-linearity parameters characterizing the bi- and 
trispectrum. Section [HI A[ Based on the formalism developed in IsoJ, we estimate /nl for narrow 
mass spectra in Section IIIIB H providing a simple approximation. Then, we take general broad 
mass spectra, properly described by the MP-distribution, to compute Jnl without any additional 
approximations. Section [HI B 2[ The resulting additional contributions to /nl due to the evolu- 



^ This formalism is a crucial extension of the linear (5 TV- formalism [39 . [sTj , and needed to deal with higher order 
correlation functions. 

^ There is a subtlety regarding initial conditions in A/'-fiation, since there is no attractor solution for quadratic 
potentials, making A/'-fiation sensitive to the initial field values - we come back to this issue later on. 
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tion outside the horizon are discussed and compared to previous analytic approximations, Section 
IIIIB 31 Finahy, we consider relaxing the slow roll condition in Section IIII C\ before we conclude in 
Section HYl 



II. TVA-FLATION AND SLOW ROLL 



We start by considering the action for M scalar fields, 

,2 /, 



S 



m: 



(1) 



which we assume to be responsible for driving an inflationary phase (see e.g. p] for a review on 
multi- field inflation). The unperturbed volume expansion rate from an initial flat hypersurface at 
t* to a final uniform density hypersurface at is given by ^ 



N{tc,Q 



Hdt. 



(2) 



where H is the Hubble parameter. 
In AA-flation 



28( 1 ■ the M ~ 1000 scalar fields that drive inflation are identified with axion fields. 



If one expands the periodic axion potentials around their minima all cross-couplings vanish 
Hence, in the vicinity of their minima, the fields have a potential of the form 

M 

W{ipi,ip2,...,w) = ^Va{(Pa) 
A=l 



El 

A=l 



(3) 
(4) 



We arranged the fields according to the magnitude of their masses, that is 



7714 > fJiB A> B. It 



IS 



2C 



21 



22] ^, where 



should be noted that AA-flation is a specific realization of assisted inflation 
the many scalar fields assist each other in driving an inflationary phase, so that no single scalar 
field needs to traverse a super-Planckian stretch in field space. 

Further, the spectrum of masses in (jl]), which were assumed to be equal in 2a], can be evaluated 



more accurately by means of random matrix theory and was found by Easther and McAllister to 
conform to the Marcenco-Pastur (MP) law 



1 



311 ]. This results in a probability for a given mass of 
0, (5) 



* Af refers to the number of scalar fields, whereas is the number of e-folds. 

^ See also [23I] for another realization of assisted inflation based on M-theory from multiple M5-branes. 
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where /3 and a completely describe the distribution: o is the average mass squared and /3 controls 
the width and shape of the spectrum (see Fig. [T]). As a consequence, the smallest and largest mass 
are given by 



a = o' 



(6) 



m^ = 6 = a2(l + v^)2. (7) 

In AA-flation, /? can be identified with the number of axions contributing to inflation divided by 
the total dimension of the moduli space (Kahler, complex structure and dilaton) in a given KKLT 



compactification of type IIB string theory 



3l[. We carry out the remainder of this analysis by 



treating /3 as a free parameter, but keeping in mind that /3 ~ 1/2 is preferred, due to constraints 



28l |. We will not need to specify cr, since it 



arising from the renormalization of Newton's constant 
will cancel out; however, its^magnitude is of course constrained by the COBE normalization just 
as rr? in chaotic inflation [25]. 

At this point, we introduce a convenient dimensionless mass parameter 

2 

= ^ , (8) 



as well as the suitable short-hand notation 



z s V/S. (9) 
= (10) 



Expectation values with respect to the MP-distribution can then be evaluated via 

1 ^ 

^2 l< 



Since /(a;) = x'^^^y^^ will appear frequently in our analysis, we introduce a more convenient 
notation and define functions namely 

Tl{y)^T^{y^'^), (13) 

where 

To.{y,)^ C ^{\-s){s-i-^)s''u,'ds, (14) 
J\li 
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(a) p{,r / < X >} 



p{xl < X >) 




0.5 1 1.5 2 



xl < X > 



FIG. 1: Probability of a given mass according to the Marcenko-Pastur distribution from ([5]), depending on 
l3 and the dimensionless square mass x = rn^ /ra\, rescaled with the expectation value (x) (also dependent 
on (3): (a) 3D-plot for 0.1 < /3 < 1, (b) shces for (3i = 1/4, (32 = 1/2 and = 3/4; the closer /3 is to one, 
the broader the mass spectrum becomes. 



so that the expectation values become 



27rz2 



(15) 



A few useful properties of the ^^-functions and analytic approximations can be found in appendix 

m 

Throughout most of the analysis, we restrict ourselves to the slow roll approximation. As 
explained above, N fields contribute to the energy density of the universe through a separable 
potential. In this regime, the dynamics of AA-flation are as follows: firstly, note that the field 
equations and Friedman equations can be written as 

dVA 



3H^ 



-V'a. 



W . 



Here and in the following we set the reduced Planck mass to nip 
mation is valid if the slow roll parameters 



ivl 

2 VF2 



w ' 



(16) 
(17) 

57rG)^^/2 = i_ This approxi- 

(18) 



are small (e^ <C 1, ?/i <C 1) and 



A/" 

e = ^ < 1 



A=l 



(19) 
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holds. The number of e-folds of inflation becomes 

/•c AT 

N{t,,Q = - J^TT^dvPA, (20) 

J* j=i A 

and the field equations can then be integrated to yield 

,.c / c \m\/ml 

^ = ( ^ ) • (21) 

Notice that this relationship between fields does not correspond to an attractor solution. As a 
result, predictions of AA-fiation can depend on initial conditions - admittedly, a less attractive 
feature of the proposal. 

There is another subtlety of AA-fiation: if the mass spectrum is broad, that is ^ 1 correspond- 
ing to the limit /3 ^ 1 , the heavier fields will drop out of slow roll early, even as infiation continues 
(this is the case for the preferred value of /3 ^ 1/2, corresponding to ^ ~ 34). The effect of these 
fields on non-Gaussianities cannot be estimated properly by using the 5A^-formalism, which we 
employ in Section IIII Al and IIIIBi We come back to this issue in Section IIII CI 



III. NON-GAUSSIANITIES 



Recent observations of e.g. the cosmic microwave background radiation (CMBR) 



or the large scale structure of the Universe 



3,3 



flflflQ 



made it possible to measure two point statistics 
of scalar perturbations to high accuracy. Existing measurements are consistent with a Gaussian 
spectrum, for which all odd correlation functions vanish, while the even ones can be expressed in 
terms of the two point function. For instance, the non-linearity parameter /atl, a measure of the 
three point function or bispectrum, is constraint to lie between — 54 < /tvl < 114 by the WMAP3 
l|. Future CMBR measurements, e.g via the Planck satellite [l^j, are expected to 



data alone 



tighten this bound up to | /nl |~ 5 (See eg. [40| for a review, and [4J| to relate a given primordial 
bispectrum to the one imprinted onto the CMBR). Furthermore, an all sky survey of the the 
21-cm background in the frequency range from 14 MHz to 40 MHz with multipoles up to 10^ 
could potentially limit this parameter down to | /atl |~ 0.01 [16]; but it should be noted that this 
proposal is highly optimistic. If such an observation is indeed possible, higher order correlation 
functions might also be within reach of observation. 

It is expected that more complicated models of the early Universe, such as A^-flation, lead to 
larger contributions than simple single field models. The physical reason for this expectation is the 
presence of isocurvature modes, which in turn cause evolution of fiuctuations even after crossing 
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of the Hubble radius (referred to as horizon crossing in the following), see e.g. 43]. Therefore, we 
will give a thorough examination of non-Gaussianities in models of AA-flation, only restricted by 
the slow roll approximation. 



A. Using Slow Roll and the Horizon Crossing Approximation 



To estimate t 



formalism 
in 



34 



le magnitude of the three and four point functions we use the non linear 6N- 



35[ , which is an extension of the linear (5A^- formalism first proposed by Starobinsky 



3a] and extended by Sasaki and Stewart [37] among others 



33 



441 ]. In this approach one 



relates the perturbation of the volume expansion rate 6N to the curvature perturbation which 
is possible if the initial hypersurface is flat and the final one is a uniform density hypersurface 3^]. 
Sfotice that this quantity is conserved on large scales in simple models, even beyond linear order 
45I, 



461 ] ^. Given this relationship between the curvature perturbation and the volume expansion 



rate, one can evaluate the momentum independent pieces of non-linearity parameters, which are 
related to higher order correlation functions, in terms of the change in during the evolution of 
the Universe, see e.g. 3]. 

Following this approach, the power-spectrum and the bispectrum i?^ can be computed. 
The ratio or to is proportional to the non-linearity parameter /nl, modulo a momentum 
dependent prefactor and factors of 2tt. This computation was performed in 
result 



33 



43 



4J] with the 



tJnl = 7^(1 + /J H ; : 72 ' 

5 16 f^M ,.2^ 



N' 



c 



16 r^JNLi 



where we introduced the short hand notation 

Na = 
Nab = 



dN 



(22) 
(23) 



(24) 
(25) 



^ The separate Universe formalism developed by Rigopoulos and Shellard in e.g. [40] is equivalent to the ^A'^- 
formalism. 
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and we refer the interested reader to [33], [43] for details. We shall focus on the second term in (j22p 
since the first one is known and rather small ^. If we further use the summation convention for 
capital indices, we arrive at ^ 



6 (4) ^ NaNeN^"" 



PnI = ,t In.. ■ (26) 



To estimate the magnitude of the four point function, we use the momentum independent 



parameters rjyi and g^L sts introduced in 



431 (see also |4g|]), 



= {NoND)S ' (27) 

25 NabcN'^N^N^ 
= 54 (iV^Ar^)3 • (2^) 

As mentioned above, we use the horizon crossing approximation in this section in order to 
compute the derivatives of the volume expansion rate; that is, we assume that modes do not evolve 
once they cross the Hubble radius at t*. With this in mind, we can set f'^ = ^. Since we also 
use the slow roll approximation, we can write ^ as 

fc ^ Y 

Nit,,Q = - / ^-^dipA. (29) 
A=l 

Confining ourselves to AA-flation, that is, if we make use of Va = employing equal 

energy initial conditions, m\Lp*^ = m?^Lp*Q , we arrive at 

^AB^ACflAyf^ZA ( 
-2 

'A 

where all potentials and slow roll parameters have to be evaluated at t*, and we replaced first 
derivatives with respect to ^\ by = ^j2eAW . 



Nabc = - "^/^ ( ""^''^ + ) ' (^2) 



^ On geometrical grounds, we know < / < 5/6 [itI. l4a| (/ characterizes the shape of the momentum triangle and 
is largest for an equilateral triangle), while r is the usual tensor:scalar ratio. The observational upper limit on this 
quantity depends on the priors used in the fitting process, but we can reliably estimate that r/16 < 0.1 This 
bound will be improved by future experiments, for example Clover 48], among others. 

* Note that /ml was already computed in [23] within the horizon crossing approximation, but we repeat it here for 
completeness. 

^ Note that some authors prefer to set ip"^ — (p\ instead - the different choices lead to the same results if modes are 
indeed frozen after horizon crossing. 



10 



It is now straightforward to evaluate the non-Gaussianity parameters to 

_6 (4) _ l_ . 
^Jnl - 2iv' ^-^^^ 

^NL = , (34) 

9NL = 0. (35) 

These parameters are leading order in the slow roll approximation, independent of the mass spec- 
trum of AA-flation and indeed indistinguishable from a single-field model with quadratic potential 
This is expected, since we neglected the main feature distinguishing multi-field models from 
single-field ones: the evolution of perturbations after horizon crossing due to the presence of isocur- 
vature modes. The vanishing oi qnl is due to the use of quadratic potentials so third derivatives 
of the potentials vanish. 



B. Beyond the Horizon Crossing Approximation 



We saw in the previous section that AA-flation is indistinguishable from single-field models with 
respect to non-Gaussianities if the horizon crossing approximation is used. Consequently, we must 
incorporate the discriminating feature ~ the evolution of modes after HC. To simplify matters, we 
make further use of the slow roll approximation, and restrict ourselves to computing f^\, which 
is the most observationally constrained parameter. 



The general expression of was computed in 



M to 



6 (4) 
^JNL 



where ua is given by 



UA 



AVa ^ e\ 



(36) 



W* W* 

with AVSi = V^ — V^>0, and the symmetric ^-matrix 



^2 



1-^ 



.A=l 



(37) 



(38) 



To evaluate this expression we first compute the field values at t^, and (see [391] for details). 
To do so, we need 2M conditions, given by the N — 1 dynamical relations between the fields from 

m 



9 / 9 

c \ mi/m . 



(39) 
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— 1 initial conditions, chosen for simplicity to be equal energy initial conditions 



10 



2 



^* = ^V\, (40) 



a condition that stems from ^ using the requirement that t^, be e-folds before t^. 

m = Y, Mf - , (41) 

1=1 

and the last one by demanding that slow roll ends for at least one field at tc- In the present case, 
the field with the largest mass leaves slow roll first 39| when 



VM = ^- (42) 

Once the masses are specified, solving these conditions is usually possible. We distinguish two 
cases: the first one involves narrow mass spectra that have (3-^1. This will result in a simple, 
analytic expression; the second one deals with more realistic broad mass spectra, e.g. for (3 ~ 1/2. 
To simplify our notation, we suppress the superscript c from here on. 



1. Narrow Mass Spectra 

By a narrow mass spectra we mean 



2 

mf 



5a = 1-^«1, (43) 
1 ^ 

^^JfY.^^-^^- (44) 

A=l 

Obviously, this case corresponds to the limit /3 ^ 0, which is not theoretically favored, but is, 
nevertheless, the easiest one to consider. Using the Marcenko-Pastur distribution, we arrive at the 
relation 

5 = 1- (x-i> (45) 

— ^> (-) 

between 5 and z = y^, where we made use of appendix [A] in the last step. 



^'^ Since there is no attractor solution for quadratic potentials, there is an unavoidable dependence on initial conditions 
in A/'-flation. 
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Using (P2]) - (|i6]) one can evaluate the field values ip\ and (p\ and the corresponding slow roll 



parameters, which can be found in 
becomes 



39( 1 ■ There, it was also shown that the second term in (j36|) 



'-t:?' = 0(*ViV=), (47) 



that is, it is second order in the slow roll parameters, as well as second order in 5. 
Using the same method, the first term in ()36p becomes 



,,2 



E..r,t(l-.l^j_ 1 f,_^\^oiS^), (48) 



2{2N + 1) V 2N + 1 



which includes, naturally, the contribution proportional to already present in the horizon 

crossing approximation. Hence we arrive at 

where we only kept the leading order contribution in 6. If we now use 5 from (j46l) and 

^^ = 77^' (50) 



1 + ^ 



we get after expanding in z = \/]3 



|/a = 7^fl-^l+0(,J). (51) 



(2Ar + i) \^ 2N + 1^ 

This is our first major result: if we compare the above with (j33p . we observe an additional 
term proportional y^, which vanishes if all the masses coincide. This is expected since, in the 
equal mass case, there is no evolution of modes after they cross the horizon. However, isocurvature 
modes get sourced for a non-zero width of the mass distribution. This in turn causes modes to 
evolve even after their wavelength becomes larger than the Hubble radius, leaving an imprint onto 
fifi- Note that the correction term in (j5ip carries an additional slow roll suppression As a 
consequence, even if experiments improve so far as to probe non-Gaussianities created in single- 
field inflationary models (which are of order 1/A^) we will not be able to distinguish them from 
AA-flation by means of their primordial non-Gaussianity, as long as non-Gaussianities are generated 
during slow roll and the mass spectrum is narrow. For simplicity, we restricted ourselves to small 
P, so that expanding in terms of this small parameter is possible; hence, the resulting correction 



See also |50|] for non-Gaussianities in an exactly soluble (toy) model of multi-component slow-roll inflation. 
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is even more suppressed by y^. Now, considering that the preferred mass spectrum in A^-flation 
corresponds to /? ~ 1/2, which is indeed broad, we must compute the exact expression for /^|^ 
in the next section. Since the resulting, cumbersome expression is vahd even for small values of 
/3, we can compare it with the straightforward analytic expression in (jSip (see Fig. [3]and|3|). It 
should be noted that the contribution due to (j47|) will also be negligible for broad mass spectra, in 
agreement with the conclusions of 3S|. 



2. Broad Mass Spectra 

In AA-flation the mass spectrum is known to conform to the Marcenko-Pastur distribution ([5]). 
Armed with the expectations values introduced in (jlip . along with (|15p we can evaluate all sums in 
()36p . However, before we can perform this replacement, we need to proceed in a similar manner as 
we saw in the previous section, where we evaluated the narrow mass spectrum. Strictly speaking, 
we have to calculate the field values at and t* from (|39p - ()42p in order to compute the potentials 
and slow roll parameters that appear in ([36]) - (j38]) . 

Before we continue, we should mention another subtlety in our analysis: we take to be the 
time when the field with the largest mass leaves slow roll ()42p : this time does not correspond to the 
end of inflation if the mass spectrum is stretched out considerably; since there is no cross coupling 
between the fields, the remaining fields can successfully continue to drive inflation even if a heavy 
field leaves slow roll. Consequently, several e-folds of inflation should be expected to follow after 
f^. Henceforth, the volume expansion rate appearing in our expressions may be smaller than the 
usual N M 60. We will discuss the effects of fields which leave slow roll while inflation continues in 
Section ImCl 

Given this caveat, let us now evaluate the fleld values ip\ and (^9^: plugging ([39]) and (00]) in 
(|4ip we arrive at 

*2 _ 4A^ 1 

- Ar(x-i)-(x-v) ^ ^ 

where we replaced the sums by expectation values as introduced in ([TT]) and defined 

..f. (53) 



Exact in the statistical sense using the large J\f limit, since we are going to use the Marcenko-Pastur mass distri- 
bution in order to evaluate expectation values. 
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y ] (a) 
0.08 

0.06 

0.04 

0.02 



Cb) 



-8 -6 



0.2 0.4 0.6 0.8 



FIG. 2: Solving ^ numerically leads to •^l/'-pl^ = y{P) for (a) -9 < log^iP) < -1, (b) 0.1 < /3 < 0.9. 
We took iV 60 in all plots. 



On the other hand, (j42p becomes 



(54) 



2N - (x-iy^) ' 

after using ([39]) and ([^0]) as well as (f52]) . Equation ([Ml) is now uncoupled and needs to be solved 
for y. Using the definitions of the ^-functions in (|13p . we can write ()54[) also as 

•^li(y) 



= 1 - 2iV- 



(55) 



-^2- -^12(2/) ■ 

Unfortunately, one cannot solve (j55|) analytically, but it is easily done with standard numerical 
routines as implemented in e.g. MAPLE. If we denote the solution to ()55p by y(/3) (see Fig. [2] for 
a plot of y over (3 for = 60), we arrive from (139l )-(l42]) at the desired field values 

1 m 



"Pa 



2 y 

"Pa 



XA 



AN 



XA aN{T^_^-TU{y)) 



where we defined 



In the following, we suppress the argument of the ^-functions, since it is always given by y. 
It is now straightforward to evaluate the slow roll parameters appearing in (jlSp to 

* ^ f -T7O -r-l \ 

Va = \^-2 - ^-2) , 

XA 

VA = 



(56) 
(57) 

(58) 



(59) 
(60) 
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and UA from (l371) to 



^1 = ^ , (61) 
2N 

EA = XAV'"'^ r n ^ n T^, (62) 

^ N^a (^°2 - ^-2) 



UA = ^{l-r^+cxAr^) , (63) 



where we defined 



After some more algebra, we can evaluate the two components of /^]^ in ([36ll to 
2^A=1 7^ ^AIT 



fip) ^ (65) 



and 



with 



Q 



4iVa 



(66) 
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2 



^"-■^ + (68) 
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+ c^e^^f ) 
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Q ^ 


2^12 + ^-2 + 2ce [^li - 
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+ c-2^2^0') 



(69) 
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(71) 
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(75) 
(76) 
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FIG. 3: — /^]^(2A^+l)6/5 over log(iQ) (/3) computed using: a. the horizon crossing approximation — /^|(2A^+ 
1)6/5 — 1, b. the (5-expansion from (jSip . c. the "exact" expression from (|78p and d. the approximation from 
((80)) . We took = 60 in all plots. Note that b. and d. are both good approximations up until (3 ^ 0.1. 

so that 

= 2 (/(/?) . (78) 
This is our second major result and will be discussed in the next section. 



3. Discussion 

A plot of f^l^ over f3 can be found in Figures [3] and U where it is also compared with the 
analytic approximation in (15ip . the horizon crossing approximation — /^|6/5 = 1/{2N + 1) and 
the approximation in (jSOp . Throughout the analysis, we took = 60 as a rough estimate for the 
number of e-folds. It is evident that the approximation in (j5ip is good up to /3 < /3 ~ 1/10. In 
this region, the leading order contribution to the exact expression in (j78p stems from the prefactor 
in ([67|) . which includes a dependence on /3 via Q defined in (j69|) . and the first summands, so that 
one may also use 

^ (80) 



2N 

as an approximation for small (3. Here = (1 ~ ■2)^/(1 ~ z'^) from (lASh and J~^2iy) defined 
in ()13p where y(/3) is the solution to (I55p . Naturally, we recover the horizon crossing result (|33p in 
the limit that /3 — > 0. 
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FIG. 4: —/nl^/Si over (3 computed using: a. the horizon crossing approximation — /^|(2A^ + 1)6/5 — 1, 
b. the (5-expansion from (fSTj) . c. the exact expression from ([78]) and d. the approximation from (l80|) . We 
took = 60 in all plots. Both approximations fail to recover the turn of observable in Figure (b). 

Both, the (5-expansion and the above approximation in (jSOp . fail for /? ~ 1/2, see Fig. [H 
However, the contribution due to F defined in (j67|) is neghgible even for very broad spectra (e.g. 

: there, two-field models with 



up to /? = 9/10 in Fig. |3]b), in agreement with the conclusions of 
a large ratio of the two masses were solved analytically and an additional slow roll suppression was 
found for F. Hence we may use 

f(4) ^ (RU 
^ 6ANa ^ ' 

(.^£2 - 2^i2 + ^-2 - {:f\ - + c?enf 

as an approximation in the region that is preferred in AA-flation. 

Around the preferred value of /3 = 1/2 the magnitude of —fj^i is smaller than the horizon 
crossing result, but only by a few percent (see Fig. Such a deviation will never be observable. 
The minimum is reached for /3 ~ 0.74 and —f^\ increases for larger values of [5 so that it catches 
up with the horizon crossing result around (5 ~ 0.88. For even larger values of [3 the magnitude 
of the non-linearity parameter increases more and more, seemingly becoming significantly large. 
However, one should take this result with caution, specifically the limit 13 ^ 1 Here, we took 
the final time tc to be the time at which the heaviest field leaves slow roll. In addition, we assumed 
that sixty e-folds of inflation occurred between i^, and tc- However, if the spectrum of masses is 



Remember that the limit /3 — > 1 corresponds to an infinitely broad mass spectrum, so /3 < 1 always. 
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indeed very broad, there will be a considerable amount of inflation even after the heaviest fields 
leave/left slow roll; hence, the potentially large value for —f^i at might very well be a transient 
phenomenon, due to a few heavy fields. Note that according to the Marcenko-Pastur distribution 
the majority of fields will have relatively light masses in the broad spectrum case, see Figure [J and 
that the majority of the masses are smaller than the average one for /? close to one. As a result, 
one might actually neglect the few heavy fields altogether, meaning, one might want to truncate 
the mass spectrum, since heavy fields will roll towards their minimum a lot faster than light ones. 

Of course, if a field leaves slow roll and starts to evolve faster, our formalism is not applicable any 
more up until the field settled in its minimum. We propose a method to estimate the production of 



5l| for recent 



non-Gaussianities analytically during these intervals in the next section (see also 
numerical work). 

Before we continue, we would like to remind the reader of the limitations of our approach: 
first, we focused on potentials without any cross coupling between the fields. One can argue 
in favor of vanishing couplings in the case of AA-flation if the fields stay close to the minima of 
3ll ]. but in general such an assumption is rather artificial If such couplings 



their potential 



are present, one should expect an enhanced production of non-Gaussianities. We also considered 
only quadratic potentials with mass spectra that conform to the Marcenko-Pastur distribution. 
We focused on this class, since the MP-distribution is expected to properly describe the spread 
of masses in AA-flation in the large M limit. We do not expect qualitative differences for other 
spectra. Nevertheless, if the potentials are not quadratic but quartic or exponential, we expect 
an additional suppression: since an attractor solution is present for potentials of these types (see 



e.g. [23]), isocurvature perturbations will be suppressed and in turn any evolution of modes after 
horizon crossing will also be suppressed, resulting in an additional reduction of non-Gaussianities. 
Lastly, we considered equal energy initial conditions, mainly to simplify the computations. Since 
there is no attractor solution for quadratic potentials, there is a dependence on the chosen initial 
state. This unavoidable sensitivity to the initial configuration of fields may be considered a flaw 
of AA-fiation, since the model becomes less predictive. However, the evident slow roll suppression 
of non-Gaussianities is insensitive to the chosen initial state, see e.g. the two-field cases studied in 



^"^ We thank F. Quevedo for useful comments regarding this point. Naturally, this means that the conditions for 
successful assisted inflation are hard to satisfy. 



19 



C. Beyond the Slow Roll Approximation 



We saw in the previous sections that multi-field inflationary models like AA-flation do not gen- 
erate large non-Gaussianities during slow roll. Henceforth, one can hardly discern them from 
simple single-field models. We expect our conclusions to be quite general, applicable to any kind 
of multi-field model during slow roll, given that the potential is separable and the kinetic terms 
of the scalar fields are canonical. However, if fields are evolving faster, one might still produce 
measurable contributions to the non-linearity parameter /^vl or higher order parameters such as 

QNL or TNL- 

Consider for instance the case where a few fields start to evolve faster while inflation still 
continues. Whenever a field behaves that way, the trajectory in field space will make a sharp turn 
and consequently, isocurvature perturbations will cause the adiabatic mode to evolve so that non- 
Gaussianities are generated. This additional source of non-Gaussianity could be recaptured by an 

effective one-field model: first, replace the multiple inflaton fields by one effective degree of freedom 

n 

o", which evolves due to an effective potential V{a) [8]. Doing so corresponds to the horizon crossing 
approximation at the perturbed level. Anyhow, this should be a good approximation since we know 
already that the additional contributions due to the evolution of modes after horizon crossing will 
hardly ever be observable. Further, assume that one of the fields leaves slow roll and starts to 
evolve faster. To incorporate the effect of the kink in the multi-field trajectory, introduce a step 
in the effective potential, followed by a subsequent re-definition of the effective model where the 
field under consideration is omitted. Comparing the effective model before and after removal of 
the field will provide the step width Act and height Al/. Hence, we arrive at a toy model with a 
single inflaton field, but sharp steps in the inflaton potential. Such steps can cause considerable 



non-Gaussianities [53|], in addition to a ringing in the power spectrum, which might actually be 
easier to observe than l itself 0, [3] • Subsequently, the fields that left slow roll will settle in a 
minimum of their potential without influencing the dynamics of the universe any more. 

Nevertheless, we do not expect this case to be realized in AA-flation: consider a broad mass 
spectrum (e.g. for (3 > 1/2) as a concrete model; here, the heavy fields violate indeed the slow 
roll condition \rii\ < 1 while inflation still goes on, but they experience an extra damping, not 
acceleration, since r] > 0. Since they evolve slower than expected from the slow roll approximation, 
they should not cause any additional non-Gaussianities, but simply hang around up until the light 
fields start to evolve faster. This will occur shortly before the onset of (p)re-heating, corresponding 
to the breakdown of slow roll for a in the effective single-field description. Consequently, we expect 



20 



53]) in TVT-flation. 



(p)re-heating, to be the primary additional source of non-Gaussianity (see e.g. 
(P)re-heating should occur in a slightly stretched-out manner during that last one or two e- 
folds, since fields do not contribute all at once. Since the universe is evolving throughout, we 
expect expansion effects to be relevant. Henceforth, simple parametric resonance models ignoring 
the aforementioned expansion will not be suitable. However, stochastic resonance [55], that is, 
broad parametric resonance in an expanding universe, seems to be a possibility. 

A careful examination of the above mentioned topics, especially (p)re-heating, is in preparation 
61|. 



IV. CONCLUSIONS 



In this article we considered AA-flation as a concrete realization of assisted inflation motivated 
by string theory. Interested in distinguishing this multi-field model from simple single-field ones, 
we computed primordial non-Gaussianities. 

After reviewing dynamics of AA-flation during slow roll, we evaluated non-linearity parameters 
characterizing the bi- and trispectrum in the horizon crossing approximation. In this limit M- 
flation, and other multi-field models, become indistinguishable from their single-field analogs. 

As a consequence, we incorporated the evolution of perturbations after horizon crossing. This 
evolution is due to the presence of isocurvature perturbations and provides a means of discrim- 
inating multi-field models from single-field ones. Focusing on the non-linearity parameter /tvl, 
which will be heavily constraint by observations in the near future, we evaluated its magnitude 
for narrow and broad mass distributions. In AA-flation, this mass spectrum is described by the 
Marcenko-Pastur law. We identified additional contributions, which turned out to be only a few 
percent of the horizon crossing result so that they are unobservable. 

The smallness of the additional terms is due to the slow roll approximation employed in this 
paper. Henceforth, we turn our attention to dynamics beyond slow roll; we argue that large, but 
possibly transient, contributions to f^L should be expected from fast rolling fields. Such fields are 
not expected in case of AA-flation, but might be present in other multi field inflationary models. 
We suggest an effective single-field model with steps in its potential to retain the main physical 
effect of these fields. 

A study of (p)re-heating in AA-flation, which should provide an additional source of non- 
See e.g. [H^, [H^, [H^, [H, [S^, [6^ for a sample on the extensive literature on (p)re-lieating. 
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Gaussianity, is in preparation 
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APPENDIX A: THE T„ FUNCTIONS 



Here we would like to gather some properties of the functions 



(Al) 



where C = (1 + ^) /{^ ~ ^) a-iid z is a free parameter °. First note that analytic expressions are 
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32| 



2{i-l) 



(A2) 



known if w = 1: for q > — 1 the functions become 

by relating T to the moments of the Marcenko-Pastur mass distribution, as analyzed in 32]. 
Furthermore, the expectation values and (a^^^) were computed in 31|], yielding 




•^-2(1) 
•^-3(1) 



1 — a ' 



(A3) 
(A4) 



(l-z2)3 a ' 
with a = (1 - 

We can also write down analytic expressions for general ui in the limit ^ — > oo, which corresponds 
to the limit z — > 1: first note that 



Tq{u}) = limj^o(^) 



Ji(ln(y)/2) 



(A5) 
(A6) 



21n(y) 

where T is a Bessel function of the first kind. All other J^q, can be computed via recursion since 



duj 
1 



^ 



(A7) 
(A8) 



It is identified with in A/'-flation. 
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follows directly from the definition ()Aip . 
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